Localized polarons and doorway vibrons in finite quantum structures 
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We consider transport through finite quantum systems such as quantum barriers, weUs, dots or 
junctions, coupled to local vibrational modes in the quantal regime. As a generic model we study 
the Holstein-Hubbard Hamiltonian with site-dependent potentials and interactions. Depending on 
the barrier height to electron-phonon coupling strength ratio and the phonon frequency we find 
distinct opposed behaviors: Vibration-mediated tunneling or intrinsic localization of (bi)polarons. 
These regimes are strongly manifested in the density correlations, mobility, and optical response 
calculated by exact numerical techniques. 



I. INTRODUCTION 

Recent progress in nanotechnology has triggered a sys- 
tematic study of electronic transport in microscopic sys- 
tems weakly coupled to external electrodesl. In such de- 
vices the active element can be a single organic molecule, 
but also a suspended Carbon nanotube, and may be 
thought of as a quantum dot contacted to metallic leads 
that act as macroscopic charge reservoirs. In small quan- 
tum dots energy level quantization becomes as important 
as electron correlations. Additionally vibrational modes 
play a central role in the electron transfer through quan- 
tum dots or molecular junctions (see, e.g., the topical 
review Ref . 12) • 

The electron-phonon (EP) interaction is found to par- 
ticularly affect the dot-lead coupling. Here electronic 
and vibrational energies can become of the same or- 
der of magnitude, e.g. when Coulomb charging is re- 
duced by screening due to the electrodes^. The same cir- 
cumstance holds in the polaron crossover regime, where 
the electrons are dressed by a phonon cloud, implying 
that phonon features for the current through the quan- 
tum device are of major importance^. Phonon and po- 
laron effects in nanoscale devices have been extensively 
discussed e.g. for (magnetic) molecular transistors^i^i^, 
quantum dots- , tunneling diodes and Aharonov-Bohm 
rings* , metal/organic/metal structures^, or Carbon nan- 
otubes^". 

In this paper we study the electronic properties of 
various EP coupled quantum systems. We consider 
one-dimensional structures, where the "quantum de- 
vice" is sandwiched between two metallic wires charac- 
terized by (tight-binding) electron hopping amplitude t, 
local Coulomb interaction C7, and EP coupling Ep (cf. 
Fig. [T]). Such systems may be described by a generalized 
Holstein-Hubbard Hamiltonian. The Holstein-Hubbard 
modeliiii^ii^iii^ is not completely realistic, of course, as it 
only includes local electron-phonon and electron-electron 
interactions as well as a coupling to (dispersionless) op- 
tical phonons. However, we are interested in fundamen- 
tal phenomena arising from the combination of electron- 
phonon interaction and "confinement" in discrete quan- 



tum structures. Besides many aspects of finite (EP cou- 
pled) quantum systems may be understood using such 
simplified effective models^i^ii^. 



II. MODEL 

Allowing for site-dependent potentials and electron- 
phonon/electron interactions, the tight-binding Holstein- 
Hubbard Hamiltonian takes the form 



H = ^A,n,,-t^(4 



R.c.)+cjoY, bib. 



(1) 



Here = A -f A^, where the potentials, A^ on site i, 
can describe a tunnel barrier, disorder, or a voltage basis. 
Since we will treat left and right leads in equilibrium, we 
choose A = throughout the sample, neglecting a bias 
between the metal leads and, in order to avoid spurious 
multi-scattering from the boundaries in a finite system, 
we take periodic boundary conditions. The parameter 
Ui {Ui = U + Ui), can be viewed as additional Hub- 
bard interaction or charging energy of, e.g., a quantum 
dot molecule. The parameter gi = [{sp + £p,i)/wo]^^^ de- 
scribes the local coupling of an electron on site i to an 
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FIG. 1: (Color online) Schematic representation of model de- 
vices described by the Hamiltonian ^ . 
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internal optical vibrational mode at the same site-'^^. Here 
{ep + £p^i) denotes the corresponding polaron binding en- 
ergy, and ujo is the frequency of the optical phonoi>i^. In 
this way the model, e.g., mimics tunneling through (sin- 
gle or double) barriers (A^ > 0), trapping of electrons, 
polarons, or bipolarons at single-impurity or double-well 
sites (Ai < 0), or transport through quantum dots with 
soft dot-lead links. 

On a translational invariant lattice (A^ = A, gi = g, 
Ui = U) the Holstein-Hubbard model can be numeri- 
cally solved by variational diagonalization in the one- 
and two-particle sectors of interest here. This holds in 
the thermodynamic limit, for the whole range of param- 
eters and any dimension (for a recent review of the Hol- 
stein (bi)polaron problem see Ref.Hi). The main result 
is a continuous cross over with increasing EP coupling 
strength, from electronic quasiparticles weakly renormal- 
ized by phonons to (small) polarons or bipolarons^. De- 
pending on the value of the adiabaticity ratio a = ujo/t, 
in one-dimensional systems, the large-to-small polaron 
cross over is determined by the more restrictive of the two 
conditions A = ep/2t > 1 (relevant for a ^ 1, adiabatic 
regime) or g^ > 1 (for a 3> 1, anti- adiabatic regime)^. 

Here we address the problems of polaron/bipolaron 
formation and phonon-assisted transport for the more 
complicated inhomogeneous barrier structures and inter- 
actions described by the above Hamiltonian. 

III. NUMERICAL RESULTS AND DISCUSSION 

In our numerical work we combine exact diagonaliza- 
tion (ED) and kernel polynomial methods2ii2^ to deter- 
mine the ground-state and spectral properties. All ener- 
gies will be measured in units of t. 

A. Single-electron case 

We first consider a single electron that tunnels through 
a single quantum barrier. The barrier height is assumed 
to considerably exceed the electron half-bandwidth. Out- 
side the barrier the electron is subjected to a rather mod- 
erate EP coupling, £p = 0.5. The chosen phonon fre- 
quency Wo = 0.4 reflects an adiabatic situation. 

Figure [2] shows the behavior of the system's kinetic 
energy 

i?fe» = "5]((cLc,,+i, + H.c.)) (2) 

as the EP coupling strength is increased at the barrier 
site. Recall that both coherent and incoherent transport 
processes contribute to Ekin- Without loss of general- 
ity we assume the barrier to be located at site 4. For 
£p,4 = the barrier is almost impermeable, consequently 
the local electron density n^^i — {ni^ + riii) near zero 
at site 4. An additional local EP interaction £p^4 renor- 
malizes the on-site adiabatic potential, i.e., it leads to a 




FIG. 2: (Color online) Kinetic energy of a single electron on 
a = 8 site ring with potential barrier A4 = 4 at site i = 4. 
The main panel gives Ekin (squares), the electron density ne,4 
(circles), and the mean phonon number nph,4, (diamonds) at 
the barrier site as functions of an additional EP coupling £p,4. 
The inset shows the variation of Ekin if the potential A4 is 
lowered keeping ep,i — (triangles up), or if ep,4 is raised 
with A4 — (triangles down). 

local polaronic level shift that softens the barrier. Note 
that the kinetic energy stays almost constant until ep^4 
exceeds a certain critical value, Sp 4. At Sp 4, the mobility 
of the electron is arrested, and the charge carrier becomes 
quasi localized at the barrier site^^ . The large number of 
bound vibrational states {uphA = (&4fe4) — 10) give rise 
to a displaced oscillator state at site 4, i.e. a new equilib- 
rium state of the lattice results, which lowers the energy. 
The jump-like transition is in striking contrast to what 
is observed if we increase only the EP coupling locally 
(without having a barrier), or if we form a quantum well 
(A4 < 0) without additional EP interaction (see inset of 
Fig- ID) ■ In these cases we found a gradual transition from 
a nearly free electron to a rather immobile particle. 

The extremely sharp polaron transition is accompanied 
by a drastic change in the optical response. The regular 
part of the optical conductivity is given by 

CTregii^) = } d(cj-CJ„), (3) 

n>0 

where j = iet Ei,^(cLci+i<T - 4+laC^a) is the current op- 
erator and \n) label the eigenstates of H with excitation 
energy Un = - Eq. 

Figure [3] shows areg{uj) and the integrated spectral 
weight 

Sreg{ljj) ^ / areg{uj')duj' , (4) 

in the vicinity of the transition, where a tiny increase 
of ep_4 (of about 8 x 10"'^, from top to bottom) sub- 
stantially changes the optical spectra. While the upper 
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FIG. 3: (Color online) Optical response for the single-barrier 
system for various ep,4. Dashed lines give the integrated 

weight Sregi^tj)- 



panel resembles the optical spectra of a large polaron 
with an absorption maximum at small frequencies and a 
rather asymmetric line shape, we found a bimodal signa- 
ture near the transition point (middle panel) and finally 
the typical (almost symmetric) small polaron absorption 
just above Sp 4 (lower panel). In this manner the system 
acts as an optical switch. 

Corresponding behavior is found if we increase the bar- 
rier (voltage bias) keeping Sp^ fixed (see Fig. [4]). Again 
the transition is "discontinuous" for small phonon fre- 
quencies, where the concept of an adiabatic energy sur- 
face holds to a good approximation. At larger phonon 
frequencies non-adiabatic effects become increasingly im- 
portant. Here the EP coupling does not work against the 
(static) barrier directly and the transition softens as in 
normal polaronic systems. Furthermore, for wq ^ 1, the 
EP coupling constant 34 is reduced (i.e., although ep^4 is 
fixed we leave the strong coupling regime). 

In Fig. d] we have included the results obtained by a 
simple approximative analytical approach to the single- 
barrier problem. Assume that 
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FIG. 4; (Color online) Kinetic energy upon increasing the 
barrier height A4 at constant EP coupling. Filled (open) 
symbols denote ED (analytical) results for N = S {N = 00). 



with 
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Of course, E{R) has to be minimized with respect to R. 
Although the kinetic energy calculated in this way ne- 
glects important contributions from multi-phonon pro- 
cesses^^-, we see that El^f^ gives a reasonable estimate for 
the critical value of A4, at least in the adiabatic regime. 
In the anti-adiabatic region, i?^f„ fails to describe the 
observed continuous cross over. This is a well-known 
shortcoming of such a kind variational approaches, which 
normally yield an abrupt polaron transition in the whole 
frequency rangoii. 



Pj = a{R) el'-J'l/« (5) 

with a{R) — tanh is the probability for finding the 
particle at site j away from the barrier site i. Then, for 
the infinite system, the ground-state energy of a polaron 
with radius R, where R = 00 corresponds to the free 
electron while R — describes a small polaron localized 
at the impurity site, is given as 

E{R)^EZ + E11^ (6) 



B. Two-electron case 

Next we investigate two electrons in a single-barrier 
structure. Now, increasing the EP coupling on a barrier 
site with strong Coulomb repulsion, we found two succes- 
sive transitions, see Fig. [SJ In the first step one electron 
becomes localized at the barrier site blocking, because of 
the large U, the second one. Raising ep_4 further, both 
particles will be trapped, forming an on-site bipolaron. 
This can be seen most clearly by monitoring the density 



FIG. 5: (Color online) Kinetic energy (squares) and mean 
electron/phonon numbers (circles/diamonds) for the case of 
two electrons in a singlet state. We assumed strong Hubbard 
interaction at the barrier {U4, = 10, A4 — 2.5). Otherwise 
U = 0, Ep = 0.5, and ujq — 0.4. The insets display the weight 
of the m-phonon state in the ground state, |cmp, for several 
characteristic ep,4. 



FIG. 6: (Color online) Two electrons in a quantum double- 
barrier system (A4,6 = 2.5), where the potential of the central 
site A5 is lowered with respect to the leads (cf. Fig.[l}. Shown 
are the kinetic energy (squares) and the on-site density cor- 
relations (triangles) for Sp = 0.5, ujo = 0.4 and U = 0, Us ~ 1 
in the upper panel, while ?7 = 9, ?75 = in the lower panel. 



correlation 



d^ = + 2(r 



)) 



(9) 



as a function of the EP coupling. Ekin and nph^4 clearly 
also show this two-step transition, being related to sig- 
nificant changes of the ground-state phonon distribu- 



tiooia^, 



see insets. The comparison of data for 



A'' = 8, 10 shows that there is almost no finite-size de- 
pendence of the results. 

Finally let us consider the double-barrier quantum dot 
structure sketched in Fig. [U with two electrons in the 
system. We plot in Fig. [6] the kinetic energy and the 
particle occupation of the barrier and embedded dot sites 
as functions of the depth of the quantum well (A5 < 
0). The upper panel describes the regime of moderate 
Coulomb interaction at the dot, with [7 = otherwise. 
Here the dot is unoccupied until its potential is lowered 
below a critical value. Then the particles initially located 
together at one of the dot-lead sites are transferred onto 
the dot. In this process they change their nature from a 
bipolaronic quasiparticle to two electrons solely (linearly 
dependent) bound by the potential well (impurity). Thus 
the ground state is a multi-phonon (few-phonon) state for 
A5 > A§ (A5 < A§). If the system has a large Coulomb 
interaction everywhere, double occupancy is prohibited 
(lower panel). Then we find initially one polaron per 
barrier (lead-dot) site and only one particle tunnels to 
the dot at A§, thereby stripping its phonons away. Note 
that the mobility is enhanced in the transition region. 

This effect is even more pronounced if we suppose that 
the EP coupling acts on the dot-lead link sites only. As 
can be seen from Fig.[7]there is large jump- like increase of 
the particle's kinetic energy if the quantum well reaches 
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FIG. 7: (Color online) Quantum-well configuration with two 
electrons. Parameters as in Fig. |S] (upper panel) but now 
Ep = 0. For comparison we show Ekin for a reference system 
without any EP coupling at the links. 



Ag . At Ag the bipolaron, located at one of the dot-lead 
link sites, dissolves and the electrons can pass over to 
the dot. Clearly Ekm decreases if we lower the potential 
of the quantum dot further but note that for A5 < A§ 
the kinetic energy is still larger than for a reference sys- 
tem without EP coupled dot-lead link sites. In this way 
the local coupling to vibrational degrees of freedom of 
the barrier opens the gate for particle transmission, i.e., 
vibronic excitations play the role of "doorway states" . 

To corroborate the importance of these quantum lat- 
tice fluctuation effects we determined the optical spectra 
below, near, and above the threshold A§. The data pre- 
sented in Fig. [5] give clear evidence for (bi) polaron hop- 
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Ag, and finally, for A5 < A§, the ground state contains 
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FIG. 8: (Color online) Optical conductivity for a "soft-linked" 
quantum dot system with two electrons (parameters as in 
Fig. [71 . The insets give the corresponding phonon distribution 
functions. 



ping transport for a shallow quantum well, with dom- 
inant phonon emission and absorption processes, but 
resonant vibration-mediated tunneling takes place for a 
deeper well. 

We emphasize that the increase of Ekm in passing be- 
low Ag is accompanied by a decrease of the total inte- 
grated weight Sregioo) of the regular (incoherent) part 
of (Treg('^) (compare dashed lines in Fig. [8] from top to 



bottom). Thus, exploiting the f-sum rule, 

- i?fcm/2 = D + Sregioo) , 



(10) 



we can conclude that the coherent contribution (Drude 
part D) to Ekin is amplified. The insets substantiate this 
interpretation. Starting from a Poisson-like distribution 



only zero- 
weight. 



one- and two-phonon states with substantial 



IV. SUMMARY 



of I 



a second maximum develops at to = for A5 > 



To conclude, investigating finite quantum structures 
coupled to vibronic degrees of freedom in the frame- 
work of a generalized Holstein-Hubbard Hamiltonian, we 
have demonstrated that interesting new physics, such 
as intrinsic (bi)polaron localization or phonon-assisted 
transmission, emerges when the energy scales set by ex- 
ternal potentials, Coulomb and electron-phonon inter- 
actions become comparable. In this regime the inter- 
play between the linear effects resulting from the barri- 
ers/cavities and the nonlinearity inherent in a discrete 
interacting electron-phonon system is of major impor- 
tance. A general understanding of vibrational effects in 
(molecular) quantum transport, however, is still far off. 
Our objects in view will be to study (i) how polaronic 
quasiparticles time evolve when passing through phonon- 
coupled nanoscale structures and (ii) how finite temper- 
ature (heating) affects the balance between coherent and 
incoherent transport mechanisms. 
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